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In order to investigate the universality of magnetic turbulence in space plasmas we analyze seven 
time periods in the free solar wind of different origin, slow or fast, and under different plasma 
conditions. The orientation of magnetic field to the flow velocity was always quasi-perpendicular. 
Unique combination of three instruments on Cluster spacecraft which operate in different frequency 
ranges give us the possibility to resolve spectra up to 300 Hz. We show that spectra measured under 
different plasma conditions have a similar shape. Such a quasi-universal spectrum consists of three 
parts: two power laws and an exponential domain. At MHD scales, Kolmogorov's law ~ k~^^^ 
is found. At scales smaller than the ion characteristic scales, a k~^'^ law is observed. At scales 
fcpe ~ (0.1 — 1), where pe is the electron gyroradius, the magnetic spectrum follows an exponential 
law exp(— A;^/^), indicating the onset of dissipation. This is the first observation of an exponential 
magnetic spectrum in space plasmas. We show that among several spatial kinetic plasma scales, the 
electron Larmor radius plays the role of a dissipation scale in space plasma turbulence. 

PACS numbers: 52.35.Ra,94.05.-a,96.60.Vg,95.30.Qd 



Space plasmas are usually in a turbulent state, and 
the solar wind is one of the closest laboratories of space 
plasma turbulence, where in-situ measurements are pos- 
sible thanks to a number of space missions 3. These 
measurements obtain time series which provide access to 
frequency spectra or to spectra of wave vectors along the 
flow. It is well established that at MHD scales (below 
~ 0.3 Hz, at 1 AU) the solar wind turbulent spectrum 
of magnetic fluctuations follows the Kolmogorov's spec- 
trum '-^ /~^/^. However, the characteristics of turbulence 
in the vicinity of the kinetic plasma scales (such as the in- 
ertial lengths Aj_e = c/wpi^ej c being the speed of light and 
LUpi^e the plasma frequencies of ions and electrons, respec- 
tively, the Larmor radii pi,e and the cyclotron frequencies 
i^ci.e = ^B/rrii^f,) are not well known experimentally and 
are a matter of debate. It was shown that at ion scales the 
turbulent spectrum has a break, and steepens to ~ 
with a spectral index s that is clearly non-universal, tak- 
ing on values in the range —4 to —2 [2, 3]. These indices 
were obtained from data that enabled a rather restricted 
range of scales above the break to be investigated, up to 
~ 3 Hz. It is not known whether such indices persist at 
higher frequencies. At electron scales, the observations 
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are difficult and our knowledge is very poor. Denskat 
Q using Helios data obtained high resolution magnetic 
spectra at 2 distances from the Sun: up to 50 Hz at 
1 AU, and up to 470 Hz at 0.3 AU. However, in both 
cases, the electron characteristic scales were not reached. 
It was only with Cluster observations that these electron 
scales were reached. For the solar wind downstream of 
the Earth's bow-shock, it was shown that the turbulence 
spectrum changes its shape around kXe ~ kpe ~ 1 [5|. 
This result was recently confirmed in the upstream solar 
wind magnetically connected to the bow-shock Q . How- 
ever, in both studies the plasma /3 (the ratio between 
plasma and magnetic pressures) was 1 and so it was 
not possible to separate the roles of pe and Ag. 

Measurements of solar wind turbulent spectra in the 
vicinity of ion and electron plasma scales may clarify our 
understanding of the processes of dissipation (or disper- 
sion) of turbulent energy in coUisionless plasmas. A num- 
ber of processes may be considered at these scales: cy- 
clotron damping at fni and fee of Alfven and whistler 
waves, respectively [7[; scattering of oblique whistler 
waves at fd < f < fee ! linear dissipation of kinetic 
Alfven waves at l/p^ < A: < 1 /pe [§, [l3| . 

In this paper we use the Cluster spacecraft data 
to analyze the free solar wind of different origin, fast 
and slow, and under different plasma conditions. While 
Sahraoui et al. [H use the FluxGate Magnetometer 
(FGM) and STAFF-Search Coil (SC) ^ at the 
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burst mode, which allow them in principle to investigate 
turbulence spectra up to 180 Hz, we complete these in- 
struments with the STAFF-Spectra Analyzer (SA), en- 
abling us to increase considerably the upper frequency 
limit, up to 4 kHz. However, as was shown in [Q], above 
100 Hz the instrument noise becomes a significant issue 
which we take into account in our analysis. 

As suggested in [ISj] , we use measurements in the mag- 
netospheric lobe (precisely, the data on 5 April 2001, 
06:00-07:00 UT) as the noise level of the instrument. The 
final spectra were obtained by substracting the lobe spec- 
trum from the solar wind spectra. A similar procedure 
has been applied by Lin [1J| for the Ulysses spacecraft 
data. The maximal frequency in our analysis is defined 
as the highest frequency where the measured spectrum is 
higher than twice the lobe spectrum, before subtraction. 

We select seven time intervals of 42 minutes when Clus- 
ter was at apogee (19 Earth radii) and spent one hour or 
more in the free solar wind: the electric field data at the 
electron plasma frequency show no evidence of magnetic 
connection to the bow shock [P. Canu, private commu- 
nication, 2009]. In Table I, the dates of the intervals 
are shown as ymmdd, and their starting times are de- 
noted by ti. Average plasma parameters for the selected 
intervals are given also in this table. Magnetic field mea- 
surements were obtained from Cluster 1. Ion moments 
(density N, velocity V and perpendicular temperature 
T±i) are measured by the CIS/HIA experiment [l^ on 
Cluster 1. The ion parallel temperatures are not prop- 
erly determined in the solar wind by the CIS instrument 
[I. Dandouras, private communication, 2009]. Electrons 
are measured by the PEACE instrument 16], mostly 
on Cluster 2. One can see from Table I that the mean 
field/flow angle, Qbv, is always larger than 60°. Other 
plasma parameters are rather variable: V varies from 
~ 360 km/s to 670 km/s, the total perpendicular plasma 
beta, (3i± + Pe± = 2fionk{Ti± + Te±)/B^, varies be- 
tween 0.7 and 3.3, the Alfven speed Va € [30, 130] km/s. 
ythi,e = \/kT±i^e/Tni^e are the ion and electron perpen- 
dicular thermal speeds, respectively; pi,e — Vthi,e/^ci,e 
are the corresponding Larmor radii. During these seven 
intervals we never observe quasi-parallel whistler waves, 
characterized by a quasi-circular right-hand polarization, 
which can be captured by STAFF-SA instrument. The 
two intervals 3 and 5 display the most intense spectra and 
are observed in the fast solar wind, a few hours down- 
stream of an interplanetary shock. 

Figure [IKtop) shows the magnetic spectrum P{f) for 
interval 5. It is calculated using the Morlet wavelet 
transform, as was done in ,17]. One can clearly recog- 
nize here two power-laws and an exponential ranges: At 
low frequencies, the spectrum is ~ /~^'^ consistent with 
Kolmogorov's law. Between fd and fx. ~ /p. (where 
/a; = V/2TTXi and /p. — V/2npi), the first break appears. 
At higher frequencies, the spectrum follows an ^ f~^'^ 
law. However, at 10 < / < 100 Hz, the spectrum is no 
longer a power-law, but follows approximatively an expo- 
nential function exp(— a(///o)°'^). At higher frequencies, 



TABLE I: Solar wind parameters for selected time periods. 



Nb 


1 


2 


3 


4 


5 


6 


7 


ymmdd 


10405 


20219 


30218 


31231 


40122 


40127 50112 


ti{\JT) 


22:36 


01:48 


00:18 


10:48 


05:03 


00:36 


02:00 


B(nT) 


7.3 


7.0 


15.5 


10.9 


15.5 


9.5 


13.6 


iV(cm~^) 


3 


29 


7 


22 


20 


8 


33 


r±«(eV) 


17 


7 


40 


10 


61 


10 


14 




36 


7 


18 


16 


28 


21 


16 


V(km/s) 


540 


365 


670 


430 


635 


430 


440 


ObvC) 


85 


65 


80 


75 


85 


80 


85 




0.4 


1.7 


0.5 


0.8 


2.0 


0.4 


1.0 


l3±e 


0.8 


1.6 


0.2 


1.2 


0.9 


0.8 


1.2 


fci (Hz) 


0.11 


0.11 


0.24 


0.17 


0.24 


0.15 


0.21 


MHz) 


205 


195 


435 


305 


435 


265 


380 


Ai(km) 


130 


40 


85 


50 


50 


80 


40 


Ae(km) 


3 


1 


2 


1 


1 


2 


1 


Pi (km) 


60 


40 


40 


30 


50 


35 


30 


Pe(km) 


2 


1 


0.5 


1 


1 


1 


0.5 


K(km/s) 


95 


30 


130 


50 


75 


70 


50 


Vthi{km/s) 


40 


25 


60 


30 


75 


30 


35 


Vih.lO'(km/s) 


2.5 


1.1 


1.8 


1.7 


2.2 


1.9 


1.7 
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FIG. 1: Top: Magnetic power spectral density for interval 
5, measured by Cluster/FGM (up to 1 Hz), STAFF-SC (up 
to 10 Hz) and STAFF-SA (/ > 8 Hz) instruments in the 
solar wind. Vertical bars indicate plasma kinetic scales, where 
f\. ^ correspond to the Doppler-shifted Ai,e and fp^ ^ to pi,e- 
Power-laws f~^'^ and f~^'^ are shown. Dashed-dotted line 
indicates exponential fit ~ exp(— o(///o)°'^), with /o = fp^ 
and the constant a ~ 9. Bottom: Compensated spectrum 
by f^'^ at low frequencies, /^'^ for middle part, and by the 
exponential for high frequency part. 
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FIG. 2: (a) Magnetic spectra for 7 time periods of 42 min- 
utes; spread of fci,e for the 7 intervals is shown (b) fe-spectra 
normalized over Pq; characteristic wave numbers, kp. = 1/pi 
etc., are shown. 



/ > , the spectrum is too close to the noise level (see 
the black solid line) to draw any firm conclusions. 

To demonstrate the above scaling laws, Fig- 
ure [Tfbottom) shows compensated energy spectra. The 
low frequency part of the spectrum was compensated by 
/^■'' (solid line), the middle range - by /^'^ (dashed) and 
the high frequency part - by exp(a(///o)°'^) (dashed- 
dotted). The combined compensated spectrum is indeed 
very flat up to . 

The spectra for the seven intervals are presented in 
Figure mja). Horizontal bars indicate the spread of fci,e 
among these seven independent observations. One can 
see that the spectra have similar shapes. Their intensity 
is, however, different. To superpose the spectra, we begin 
by applying Taylor's hypothesis, which should be valid 
for the whole frequency range, as far as quasi-parallel 
whistler waves are not observed during selected inter- 
vals (as mentioned above). Thus, we assume that the 
frequency-spectra are indeed Doppler shifted /c-spectra 
P{k) = P{f )V/2n with k = 27r//y. Then, we deter- 
mine a relative intensity of the j-th spectrum, Sj with 
j = 1, ...,7, as Po{j) = {Sj/Si), where is a reference 
spectrum, and (•) indicates a mean over the range of wave 
vectors 10~^ < k < 10~^ km"-'^. With this normalization, 
the rescaled spectra may be nearly superposed as shown 
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FIG. 3; Relative spectral intensity Pq as a function of (a) 
magnetic and (b) kinetic energies; (c) Pq as a function of the 
ion cyclotron period and (d) the electron gyroradius. Linear 
fits with corresponding slopes are shown by solid lines. 



in Figure IDJb). 

One expects that the spectral level, Pq, depends on 
the solar wind kinetic, thermal or magnetic energy. The 
scatter plots shown in Figure[Ha) and (b) indicate clearly 
a power law dependance of Pq on the magnetic energy, 
and a less clear dependance on the kinetic energy, and 
thermal energy (not shown). 

To understand the meaning of the observed de- 
pendence on the magnetic energy, one may use a 
Kolmogorov-like phenomenology. Suppose first that the 
solar wind magnetic turbulence dissipates through an ef- 
fective diffusion mechanism of rjAB (rj being a - prob- 
ably turbulent - magnetic diffusivity) , and second, that 
the observed turbulence is quasi stationary. In such a 
case, there is a balance between the energy input from 
nonlinear interactions at large scales and the energy drain 
from the dissipation at small scales. This implies that the 
energy transfer rate e depends on the dissipation scale 



The depen- 



£d as e = 773£-4; thus Pq ~ £2/3 _ ^"8/3^ 

deuces observed in Figure [3{c) and (d) , Pq 
and Pq Pe^'^, are very close to the prediction of this 
phenomenological model. JVIore statistics are needed to 
confirm the observed exponents. We can state, however, 
that the observed dependences imply that pe and/or fd 
and/or fee play an important role in the dissipation pro- 
cesses in collisionless plasmas. Let us now confirm these 
results. 

From the balance between the energy input and the 
dissipation, for the Kolmogorov's spectrum E(k), it fol- 
lows as well that E{k)tdlff' is a universal function of 
Md El, El- Figure H] tests which of the kinetic scales 
is to be used as ^d to recover a universal function from 
the observed spectra. We assume for simplicity that 
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FIG. 4: Universal Kolmogorov function oc £dE{k) for hypoth- 
esized dissipation scales £d as a function of (a) kpi, (b) k\i, 
(c) kpe and (d) ///ce- 



is constant, despite the varying plasma conditions. One 
can see that the pi and Ai-normaHzations are not effi- 
cient to collapse the spectra together. Normalization on 
Ae gives the same result as for A^. At the same time, the 
normalization on pe and fee bring the spectra close to 
each other. This confirms that the electron gyroradius 
Pe and/or cyclotron periods of the particles are impor- 
tant in the dissipation. 

With the present observations it is not possible to dis- 
tinguish between pe and cyclotron periods as far as there 
is a clear correlation between pe and B. We can ar- 
gue, however, that if the cyclotron period had been the 
only dissipation scale, the turbulent cascade would have 
stopped by the cyclotron damping of Alfven waves at fd 
showing an exponential cut-off at this scale ^2d\. Solar 



wind observations show the contrary: the turbulent spec- 
trum continues up to electron scales. Thus, we conclude 
that Pe is the dissipation scale of magnetic turbulence in 
the solar wind, but we can not exclude that at fd and 
fee there is a partial dissipation by cyclotron damping. 

In the present letter we analyzed high resolution mag- 
netic spectra from MHD to electron scales. We show here 
for the first time that whatever the plasma conditions and 
the solar wind regime, slow or fast, the magnetic spectra 
have similar shape. This indicates a certain universal- 
ity, at least for the quasi-perpendicular configuration be- 
tween B and V. Such a quasi-universal spectrum consists 
of three parts: two power-laws and an exponential do- 
main. At MHD scales it follows a Kolmogorov's ^ 
spectrum, in agreement with previous observations. Be- 
tween fei and Doppler shifted Ai and pi a spectral break is 
observed. Above the break, it follows a power-law. 
At scales, of the order of kpe ~ (0.1 — 1) it follows an ex- 
ponential exp(— a(///o)^/^). This is the first observation 
of an exponential magnetic spectrum in space plasmas. 
Such spectra were predicted by the anisotropic dissipa- 
tion model of Gogoberidze Q . The author suggests that 
small scale fluctuations with oblique k are diffused on 
oblique fluctuations from the inertial range. This diffu- 
sion is anisotropic and it gives an ~ exp(— fc'^"/^) spec- 
trum in the dissipation range, where Aa = a±_ — a|| is 
the difference between the energy diffusion scaling per- 
pendicular and parallel to B. 

It is a long standing problem to distinguish between the 
role of different kinetic scales in space plasmas. We show 
for the first time that the role of dissipation scale in space 
plasma turbulence is played by the electr on gy roradius, 
as assumed by several previous authors [1, [UHH . 
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